On the toroidal nature of the low-energy El mode 
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The nature of El low-energy strength (LES), often denoted as a "pygmy dipole resonance", 
is analyzed within the random-phase-approximation (RPA) in '^"^Fb using Skyrme forces in fully 
self-consistent manner. A first overview is given by the strength functions for the dipole, com- 
pressional and toroidal operators. More detailed insight is gained by averaged transition densities 
and currents where the latter provide very illustrative flow pattern. The analysis reveals clear 
isoscalar toroidal flow in the low-energy bin 6.0-8.8 MeV of the LES and a mixed isoscalar/isovector 
toroidal/compression flow in the higher bin 8.8-10.5 MeV. Thus the modes covered by LES embrace 
both vortical and irrotational motion. The simple collective picture of the LES as a "pygmy" mode 
(oscillations of the neutron excess against the nuclear core) is not conflrmed. 



During the last decade we observe an increasing inter- 
est in low-energy El strength (LES), for a recent review 
see This interest is caused by a possible relation of 
LES to the neutron skin in nuclei and density depen- 
dence of the nuclear symmetry energy. This in turn may 
be important for building the isospin-dependent part of 
the nuclear equation of state and various astrophysical 
applications [2[ . Several different views of the LES origin 
come together. Most often the LES is interpreted as a 
" pygmy dipole resonance" (PDR) modeled as the oscilla- 
tion of the neutron excess against the nuclear core [ll, 

There are, however, serious objections against such a 
simplistic collective picture @, 0]. In fact, the landscape 
of LES may be much richer. It can embrace also the 
toroidal resonance (TR) & Q and anisotropic compres- 
sion resonance (CR) [lO, llll which both are of actual in- 
terest [l|. After exclusion of nuclear center-of-mass (cm.) 
motion, the TR and CR dominate in the isoscalar (T=0) 
channel and constitute the low- and high-energy branches 
of the isoscalar giant dipole resonance (ISGDR). Follow- 
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FIG. 1: Schematic velocity fields for the El pygmy (a), 
toroidal (b), and high-energy compressional (c) flows. The 
driving fleld is directed along z-axis. The arrows indicate 
only directions of the flows but not their strength. In the 
plot (c), the compression {+) and decompression (-) regions, 
characterized by increased and decreased density, are marked. 



ing recent microscopic studies [l2|, |l3[ , the TR dominates 
in the LES region and the CR, being strongly coupled to 
TR, also significantly contributes there. The basic flow 
patterns of these three modes are shown schematically in 
Fig. 1. The panels illustrate the PDR oscillations of the 
neutron skin against the core (a), the typical vortices of 
the TR (b), and the dipole-compressional pattern of the 
CR (c). The latter can be viewed as oscillation of sur- 
face against core and thus shares some similarity with the 
PDR picture of panel (a). Unlike the irrotational PDR 
and CR, the TR is purely vortical in the hydrodynami- 
cal (HD) sense [l^ . [isj . Thus we see that the LES can 
involve quite different flows, vortical and irrotational. 

Despite a great number of publications on PDR, TR 
and CR their possible interplay in LES was only 
occasionally discussed. There is a study within the 
quasiparticle-phonon model which discusses PDR and 
TR side by side 141. It was shown that the vortical 
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strength [15| is peaked in the LES region and the isovec- 
tor LES velocity field is mainly toroidal. Nevertheless, 
because of the dominant contribution of the neutron skin 
to the surface transition density and thus to B(E1,T=1), 
the familiar PDR picture of LES was maintained [14 1. 
Similar arguments in favor of the PDR treatment were 
earlier presented in relativistic [3| and Skyrme Q mean- 
field calculations, and taken up in most of subsequent 
publications. Recent explorations question the simple 
PDR-type collectivity of LES [1, 0, [H, though without 
analysis of LES velocity fields. 

It is the aim of this paper to give a more thorough 
exploration of the interplay and structure of low lying 
dipole modes. Not only strength functions and transi- 
tion densities but also the mode fiow patterns will be 
considered. As we will see, the actual LES flow is pre- 
dominantly of mixed TR/CR character. This conclusion 
may have far-reaching consequences for the information 
content of LES. If vorticity dominates, then only a mi- 
nor irrotational fraction of LES is relevant for the nu- 
clear symmetry energy and related problems (as was also 
worked out recently using the correlation analysis 0). 

Our study is performed for ^*'*Pb using the Skyrme 
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RPA approach with the techniques of [1^. The method 
is fully self-consistent as both the mean field and resid- 
ual interaction are derived from the Skyrme functional 
The RPA residual interaction takes into ac- 
count all the terms of the Skyrme functional including 
the Coulomb (direct and exchange) energy. The center- 
of-mass correction (c.m.c.) is implemented for the ground 
state and T=:0 dipole excitations. The parameterization 
SLy6 [2^ is used which provides a satisfactory descrip- 
tion of E1(T=1) strength in heavy nuclei [2l|. The calcu- 
lations are done in a ID spherical coordinate-space grid 
with mesh size 0.3 fm and a calculation box of 21 fm. A 
large configuration space including Iph states up to ^ 35 
MeV and additional fluid dynamical basis modes is used. 
The later allows to include global polarization effects up 
to ~200 MeV [1^, correctly extract the cm. motion, 
and fully exhaust the Thomas-Rciche-Kuhn sum rule. 

The excitation modes are first characterized by their 
strength function 

So.iEl;Lu) = 3^c^^|(*,|A/„(i;iO)|*o>|'C(^,^.) (1) 



where C(w,Wi/) = A(a;i,)/[27r[(a; - cj^)^ -I- A(cj^)^/4]] is 
a Lorentzian weight with energy- dependent smoothing 
width A(u;J_= max{0.4 MeV, (w^, - 8 MeV)/3)}, for 
details see [22|. Further, ^I^o is the RPA ground state 
(g.s.) while i> runs over the RPA spectrum with eigen- 
frcquencics uj^ and cigcn-statcs The Ma{Elfi) is 

the transition operator of the type a — {El, tor, com}. 

For E1(T=1) transitions (a =E1), we consider the 
ordinary El operator (cx rYi^^) with effective charges 
el = N/A and e'^ = -Z/A and the strength ^ is 
weighted by the energy, i.e. 1 = 1. For a= tor, com, 
we implement eg = eg = 1 for T=:0 (e^ = — e" = 1 for 
T=l) and no energy weight {I = 0). 

The TR and CR operators used in ^ read [Hi 



(2) 
(3) 



where jc{^) is the operator of the convection nuclear cur- 
rent, Yii^{r) and Yip(r) are vector and ordinary spher- 
ical harmonics. The terms with the g.s. squared radius 
(r^)o = / (Pr r'^po{r)/A account for the c.m.c, po{f) 
is the g.s. density. Note that we describe CR and TR 
operators on the same footing usin g th e current oper- 
ator. There is the direct relation [ijl Mcom(£'l/i) = 
~M'^^^{Elp)Lo / c between the CR current-dependent op- 
erator Q and its familiar density-dependent counterpart 
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FIG. 2: Strength functions calculated within RPA with the 
force SLy6. (a) E1(T=1) giant resonance The line and ar- 
row indicate the experimental width and energy centroid of 
the resonance [S^l . The pygmy region is marked, (b) Toroidal 
(solid line) and compression (dotted line) E1(T=0) strengths. 
The widths and energy centroids of the low- and high -energy 
branches of ISGDR observed in {a, a') reaction [2J| are de- 
noted, c) The same as in the plot (b) but in T=l channel. 



where p{r) is the density operator. 

The operators ©-([l]) are derived as second-order ^ 
I'^Yiij, terms in the low-momentum expansion of the ordi- 
nary El transition operator [l,[^,[i2|. Despite its second- 
order origin, TR and CR dominate in the E1(T=0) chan- 
nel where the leading cm. motion driven by the operator 
rYi^ is removed as being the spurious mode. Following 

([2|)-(l3]), TR and CR deliver information on the curl V x 

and divergence V • jc of the nuclear current. As shown in 
p^ . the corresponding velocity operators indicate that 
TR is purely vortical and CR is irrotational. 

The strength functions ([T]) are shown in Fig. 2. In 
panel a), we see a good agreement of the computed gi- 
ant dipole resonance (GDR) with the experiment [23| . 
which confirms the accuracy of our description. For the 
LES region 6-10.5 MeV (marked as pygmy), we get two 
peaks at 7.5 and 10.3 MeV in accordance with previous 
RMF calculations Panels b) and c) show TR and 
CR strengths in T=0 and T=l channels. For T=0, the 
TR and CR arc believed to constitute the low- and high- 
energy parts of the ISGDR Our results somewhat 
deviate from the experimental (a, a') data p^ ]. but a 
similar discrepancy takes place in almost all theoretical 
studies [l|. Perhaps, at 12-14 MeV not the TR but the 
low-energy CR bump (see panel b)) is observed. The dis- 
crepancy for the high-energy CR may be caused by its 
sensibility to the calculation scheme. What is important 
for our aims, the LES region should certainly host the 
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dominant and strongly peaked part of TR(T=0), the left 
flank of TR(T=1), and a non- negligible low-energy frae- 
tion of CR. In other words, we expect here a complicated 
interplay of several modes. 

To understand the LES structure, we need more de- 
tailed observables than the strength distribution. In 
the following, we will consider transition densities (TD) 
dpi,(r) — (4'i/|/5|^o) ^-nd current transition densities 

(CTD) Sj^{r) = (^'^.iJd^'o) (analogous to velocity fields). 
As we have in ^"^Pb a high density of states, it is not 
worth to look at the pattern of individual states v, which 
can vary from state to state and easily hide common fea- 
tures of the flow. Thus we will consider transition densi- 
ties and velocity fields averaged over given energy inter- 
vals. Incoherent averaging requires expressions which are 
bi- linear in the excited states I'I'i/). This is achieved by 
summing TD and CTD weighted by the matrix elements 
Dt^ = {v\DTiEl)\0) of a probe operator Dt{EI): 



Spf\r) = 
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The sums in ([S])-® involve all the RPA states \v) in the 
energy interval [wi,aj2]- Since states contribute twice 
(to D^^, and transition densities 5 pi/ 5j1), the expression 
is independent of the phase of each state as it should 
be. In dS])-®, the index /3 = p,n,0,l defines the type of 
TD or CTD (neutron, proton, T=0, T=l) by the proper 
choice of the effective charges: ~ l,e" = 0; = 



0, 



1; 



1; 



= N/A, el = -Z/A. We use 
two different dipole probe operators: the isovector Di ~ 
{N/A) J2f{rYi)i~{Z/A) J2f (rYi), relevant for reactions 
with photons and electrons, and isoscalar compressional 
Dq — relevant for {a, a') reaction. Due to 

D^^ weights, the contribution of RPA states with a large 
Dt strength is enhanced. 

In Fig. 3, the TD summed over two bins of the LES, 
6.0-8.8 MeV and 8.8-10.5 MeV, are shown. One sees that, 
up to a scale factor, the TD for the probes Di and Dq are 
rather similar, especially at 6.0-8.8 MeV. Panels a) and 
c) show that at 6.0-8.8 MeV the protons and neutrons 
oscillate in phase in the interior area 4 - 7 fm (isoscalar 
flow of the core) but neutrons dominate at larger dis- 
tances r > 7 fm (contribution of the neutron excess). 
Due to r^-factor, mainly the surface area r > 7 fm con- 
tributes to the B(E1) oc / drr^Sp{r). This would favor 
the simple PDR picture of neutron-core oscillations. At 
the higher energy, 8.8-10.5 MeV, we see mainly isovector 
motion at 6-8 fm and again dominance of neutrons at 
r > 8 fm. So, LES here is more isovector and does not 
support the PDR picture. Furthermore, the isoscalar Dq 
leads to much weaker TD at 8.8-10.5 MeV relative to the 
region 6.0-8.8 MeV. This is probably caused by the fact 
that LES at 6.0-8.8 MeV is mainly isoscalar thus gaining 
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FIG. 3: Summed r -weighted proton and neutron TD Spp^n 
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(a-b) and (c-d) at the energy intervals 6.0-8.8 MeV 

(left) and 8.8-10.5 MeV (right). 
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FIG. 4: Summed CTD (a) ijf'-'^ for GDR in the bin 11-15 
MeV and (b) fjf°^ for CR(T=0) in the bin 22-30 MeV (b). 



more from the isoscalar weight. Note that the r^-factor 
amplifies the pattern in the nuclear surface and damps 
it in the interior area at r < 4 fm. The TD in Fig. 3 
indicate that there are sizable effects in the interior. This 
will be corroborated by the flow pictures below. 

A thorough analysis of excitations should also look at 
CTD which reveal even more details than mere TD. The 
CTD for dipole states \p = 10 are presented in Figs. 
4-7. In Fig. 4, the fields for the isovector GDR (10.5- 
15 MeV) and high-energy isoscalar CR (22-30 MeV) are 
given as reference examples. They show typical GDR 
and CR fiows, see [3, [la for a comparison. In the 
CR case, the compression and decompression zones along 
the z-axis are visible, as in Fig. Ic). These plots for well 
known modes serve as benchmark and assert the validity 
of our prescription. 

In Figs. 5-7 the fields for the two LES bins, 6.0-8.8 and 
8.8-10.5 MeV, are depicted. Since D\ and fields look, 
up to a scale factor, rather similar (especially for the bin 
6.0-8.8 MeV), we will show further on only weighted 
CTD. In every figure, the actual CTD scale (common for 
all the plots) is adjusted for better view. Arbitrary units 
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FIG. 5: Proton (a), neutron (b), T=0 (c) and T=l (d) 
summed CTD in the bin 6.0-8.8 MeV. 



FIG. 7: The r -weighted proton and neutron 
bins 6.0-8.8 MeV (a)-(b) and 8.8-10.5 MeV (c)-(d). 
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FIG. 6: The same as in Fig. 5 but for the bin 8.8-10.5 MeV. 



for CTD arc used. 

Fig. 5 shows the fields in the bin 6.0-8.8 MeV. The neu- 
tron flow dominates in both interior and surface. Since 
protons and neutrons move in phase, the total flow is 
essentially isoscalar, see panels c) and d) in comparison. 
The T=0 character of the lower bin of LES is in accor- 
dance with previous theoretical results [1, [l^ and exper- 
imental flndings, e.g. for ^^^Sn [g^. What is important 
for our aims, Fig. 5 clearly demonstrates the overwhelm- 



ing toroidal flow in neutron and T=0 cases (less in the 
proton case). This is in accordance with the TR(T=0) 
strength from Fig. 2b, which is strictly peaked just at 
7-8 MeV and dominates over the CR(T=0). Therefore, 
LES at 6.0-8.8 MeV is of almost pure toroidal (vortical) 
nature! The irrotational FDR flow is not seen at all. 

The LES flelds in the bin 8.8-10.5 MeV in Fig. 6 are 
more complicated. The flow is mainly isovector in the 
nuclear interior and isoscalar at the surface, thus demon- 
strating an isospin-mixed character (again in accordance 
to RMF findings [1, [13)- Furthermore, the TR flow is 
faint. Actually there are hints of several flows: TR (b-c), 
CR (b), and familiar linear dipole (a,d). This complex 
picture reflects the fact that, following Figs. 2b and c, 
the region 8.8-10.5 MeV hosts various modes and feels 
already the vicinity to the GDR. 

Finally, Fig. 7 exhibits the r^-weighted CTD to high- 
light the role of surface nucleons (e.g. the neutron excess) 
in the peripheral reactions like {a, a') and, to a lesser ex- 
tent, photo-absorption. Flows in Fig. 7 correspond to 
the TD in Fig. 3 (a-b). The r^-weighted presentation 
weakens the interior flow and emphasizes the role of the 
neutron excess, see Fig. 7 b). Nevertheless, the LES still 
keeps their TR and mixed (TR/CR) nature in both bins, 
6.0-8.8 and 8.8-10.5 MeV. Again we cannot flnd any siz- 
able evidence for FDR flow. 

The question remains how to observe the velocity flelds 
experimentally and thus disclose the true nature of LES. 
The typical reactions mentioned above are mainly sensi- 
tive to the nuclear surface and lose the important infor- 
mation on the nuclear interior. This is especially the 
case for the isoscalar (a, a') whose response is driven 
by the operator r^Yi with a huge surface enhanced fac- 
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tor. (Note also that the most relevant {a, a') measure- 
ments of ISGDR in ^°^Pb consider the energy inter- 
val (jj >8 MeV and so, following our calculations, lose the 
TR(T=0) peaked at 7-8 MeV). Perhaps, the (e, e') reac- 
tion which can cover both nuclear surface and interior is 
the most promising tool to examine LES flows. 

In this study wc do not take into account the coupling 
with complex configurations which may be essential for 
LES [H, tl^. However the TR/CR signatures in LES 
look strong enough to be appreciably smeared out by 
this effect. 

In conclusion, Skyrme-RPA calculations have been 
performed to inspect the nature of the El low-energy 
strength (LES), often denoted as the pygmy dipole res- 
onance (PDR) and associated with the picture that the 
neutron skin oscillates against the nuclear core. Strength 
functions, averaged transition densities and averaged cur- 
rent fields (collecting contributions of all RPA states in 
a given energy interval) were used for the analysis. The 
current fields turned out to be most important to illus- 
trate the LES flows. The results show that, in agreement 
with previous studies [1, [131 , LES may be divided into 
two energy regions, 6.0-8.8 MeV and 8.8-10.5 MeV in 
our case, where the lower one is basically isoscalar and 
higher one is isospin-mixed. 



What is most interesting, LES at 6.0-8.8 McV shows a 
clear toroidal (vortical) nature while the interval 8.8-10.5 
MeV gives a mixed toroidal/compression/linear flow. No 
convincing indicator of PDR-like flow is found. This 
means that the familiar treatment of LES as the out-of- 
phase motion of the neutron excess against the nuclear 
core (arising from the analogy with light halo nuclei and 
suggested from 7'^-wcighted transition densities) is mis- 
leading. Our study does not deny the important contri- 
bution of the neutron excess to various (basically periph- 
eral) reactions. At the same time, we find that LES flow 
pattern is far from a simple PDR picture and actually in- 
volve various types of motion, irrotational (compression) 
and vortical (toroidal). In particular, LES at 6.0-8.8 MeV 
constitutes an almost pure toroidal T=0 resonance. This 
conclusion may have far-reaching consequences for fur- 
ther exploration of LES and related obscrvables. 
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